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We investigate the low temperature behaviour of the integrable ID two-component spinor Bose 
gas using the thermodynamic Bethe ansatz. We find that for strong coupling the characteristics 
of the thermodynamics at low temperatures are quantitatively affected by the spin ferromagnetic 
states, which are described by an effective ferromagnetic Heisenberg chain. The free energy, specific 
heat, susceptibility and local pair correlation function are calculated for various physical regimes in 
terms of temperature and interaction strength. These thermodynamic properties reveal spin effects 
which are significantly different than those of the spinless Bose gas. The zero-field susceptibility 
for finite strong repulsion exceeds that of a free spin paramagnet. The critical exponents of the 
specific heat Cv ~ T^^^ and the susceptibility x ~ ^re indicative of the ferromagnetic signature 
of the two-component spinor Bose gas. Our analytic results are consistent with general arguments 
by Eisenberg and Lieb for polarized spinor bosons. 
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I. INTRODUCTION 



Experiments with ultracold quantum gases are opening up exciting new possibilities for testing and exploring 
quantum effects in many-body systems (for recent reviews, see Refs. These include experiments on effectively 

one-dimensional (ID) quantum Bose gases of ®''Rb atoms in which the interaction strength between atoms is tunable 
0) S B S @]- T^^G experiments provide a striking example of realizing an integrable quantum many-body problem. 
They demonstrate the explicit fermionization of bosons and provide a direct test of theoretical results obtained for the 
integrable ID interacting (spinless) Bose gas [1, Another frontier of activity involves spinor Bose gases of alkali 
atoms in which hyperfine states comprise the pseudospins llll . Il2l ll3l . In these systems quantum collisional effects 
can produce spatio-temporal spin oscillations (spin waves) |l4l.ll5l.[l6l.[l7j. The observation of collective dynamics of 
spin waves and spin-state segregation in trapped spinor Bose gases has stimulated a wide range of interest in studying 
magnetism, topological spin defects and novel quantum phase transitions in spinor Bose gases [H, [l^ . 

Two-component spinor Bose gases have been experimentally created in a ma gne tic trap by rotating two hyperfine 
states so that the two atomic hyperfine states make up a pseudo-spin doublet [20|, e.g., the \F = 2,mp — —1) and 
\F = l,mF = 1) hyperfine states of ®^Rb. In general, spin-independent s-wave scattering dominates interactions 
in alkali atomic gases. In ID, the two-component Bose gas with spin-independent s-wave scattering can be exactly 
solved, like the spinless model, by means of the Bethe ansatz [2ll.[23|. In contrast to Fermi gases, ferromagnetic order 
emerges in spinor Bose gases as long as the interaction is fully spin independent [isl . [23j . The low-energy excitations 
of the model split into collective excitations carrying charge and collective excitations carrying spin. The charge 
excitations are phonons whereas the spin excitations have quadratic dispersion connected to spin wave excitations 
[2^ [25}. Spin dynamics in the ID ferromagnetic Bose gas have been studied recently ^2d\. Girardeau's Fermi-Bose 
mapping has been used to study the ID spinor Bose gases [l^]. In general the two-component interacting Bose gas 
provides a tunable testing ground for observing the phenomenon of spin-charge separation [2^ . 

Quantum gases with multi-spin states are expected to exhibit even richer quantum effects than their single com- 
ponent counterparts [U, 0, S| ■ Universal features appearing in the low temperature behaviour of strongly interacting 
spinor Bose gases should differ significantly from those of spinless Bose gases and the antiferromagnetic behaviour 
of Fermi gases due to their fundamentally different statistical signatures. One way to calculate the thermodynamics 
of integrable many-body systems is via the thermodynamic Bethe ansatz (TBA) [2^, [s^, [Ml, [H, [sl] , introduced by 
Yang and Yang [34] for the ID Bose gas. However, it is a challenging problem to derive exact TBA results for the 
thermodynamics of ID quantum many-body systems. Our aim here is to obtain universal characteristics of ferro- 
magnetic behaviour for the ID two-component strongly interacting Bose gas of ultracold atoms via the TBA method. 
We will see that the ferromagnetic phase associated with the spin degrees of freedom may separate from the gas 
phase in the strongly repulsive regime due to spin-charge separation. The low temperature behaviour is dominated 
by the spin ferromagnetic states, which are described by an effective ferromagnetic Heisenberg chain. In this way 
we make contact with the known results for the thermodynamics of the ferromagnetic Heisenberg chain [ssl [36j to 
derive analytic expressions for the free energy, specific heat, susceptibility and local pair correlation function for the 
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strongly interacting two-component Bose gas in terms of temperature and interaction strength. These thermodynamic 
properties reveal some novel spin effects. Our explicit results are consistent with general arguments by Eisenberg and 
Lieb [13] for polarized spinor bosons. 

This paper is set out as follows. In section [IT] we present the Bethe ansatz solution of the ID two-component 
interacting Bose gas. The ground state properties are also calculated. In section [1111 we introduce the TBA for the 
spinor Bose gas in order to study the thermodynamics at low temperatures, including the analysis of spin charge 
separation. The ferromagnetic ground state is studied by means of the solution of the TBA equations in section [TVl 
We discuss low temperature ferromagnetic behaviour for the ID strongly interacting Bose gas of atoms in section [V] 
The local pair correlation function is studied at low temperatures in section IVII Section IVIII is devoted to a brief 
summary and concluding remarks. 

II. THE TWO-COMPONENT SPINOR BOSE GAS 

The Hamiltonian describing a (5-function interacting gas of N bosons of mass m constrained by periodic boundary 
conditions to a line of length L with internal degrees of freedom is 

^ = "^IIa^+5iD S{x^-Xj). (1) 

4=1 * i<i<j<N 

For the ultracold atomic gases [37j . the coupling constant giD can be written in terms of the scattering strength 
c — 2/aiD as giD = ti^c/m. The effective ID scattering length oid can be related to the 3D scattering length 
for bosons or fermions confined in a ID geometry. The dimensionless coupling constant 7 = c/n ~ mgn)/ {h^n) is 
convenient for physical analysis. Here n = N/L is the linear density. We take 2m — h — I for simplicity in the 
following equations. However, we reinstate them where appropriate in discussing the thermodynamics of the model. 
The wavefunctions of Hamiltonian ([T]) for the spinor Bose gas are symmetric under exchange of spatial and internal 
spin coordinates between two particles. We shall see that this statistical signature triggers rather novel ferromagnetic 
behaviour in the degenerate quantum spinor Bose gas. The interaction is attractive for giD < and repulsive for 
5iD > 0. However, one should note that there is no thermodynamic limit for the attractive case in Bose gases. 

For M spin-down bosons, the Bethe ansatz equations (BAE) for the two-component Bose gas are of the form (2ll[22j 



2-lC 



N . , 1. M , . 

Xa- Kp~ ttIC T-r An, - a 



f\ Xc.-h + iic J-j-^ A„ - A^ + ic' ' 

in terms of which the energy eigenspectrum is given hy E = X^jLi ■ the thermodynamic limit, i.e., N, L ^ 00 
with N/L finite, these equations can be written as coupled integral equations in terms of the particle and hole root 
densities p(fc) and /9^(fc) (o'(A) and o-^(A)) for the charge (spin) degrees of freedom, respectively. These are 

n^ hn^ 1 1 2cp{k') „, 1 ca{X) 



27r 2t: J _Q + {k ~ k'Y 271 J _b + {k - Xf 

hr.-. 1 cp(fc) „ 1 2ca{X) 

The integration limits Q and B are determined by N/L ~ J^^ p{k)dk and M/ L = a{X)dX. At zero temperature, 

the ground state corresponds to the configuration cr(A) = p^{k) = leading to a ferromagnetic ground state. For the 
ground state there are therefore no holes in the charge degrees of freedom and no quasiparticles in the spin degrees of 
freedom. However, as the temperature increases spin strings become involved in the thermal equilibrium states. We 
shall investigate this ground state configuration via analysis of the string solutions to the TBA. 

III. THE THERMODYNAMIC BETHE ANSATZ 

For finite temperatures each of the N quasimomenta ki are real due to the repulsive interaction. However, the spin 
quasimomenta form complex strings of the form [38] A^ ^ = A" -|- i(n -I- 1 — 2j)c/2 for j = 1, . . . ,n. Here the number 
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of strings a = 1, . . . , A^„. on the real axis denotes the position of the centre of a length-n string. The number of 
n-strings satisfies the relation M = '^^nNn- It is assumed that the distribution of Bethe roots along the real 
axis is dense enough to pass to the continuum limit. After performing a standard calculation with the string solutions 
and introducing the convolution integral (/ * g){\) = /(A — X')g{X')dX' , the BAE ^ become 

OO 

CT„(A)+a^(A) = a„*p(A)- ^r„™*a™(A) (4) 

m— 1 

where 

T (W - f «|n-m|(A) + 2a|„_™i+2(A) + . . . + 2a„+TO_2(A) + a„+„(A), for n ^ m, , , 

inm[A) - I 2a2(A) + 2a4(A) + . . . + 2a2„-2(A) + a2„(A), for n = m, 

and 

1 nc 
°"(^)^2^(nc/2)2 + A2- 

The equilibrium states at finite temperature T are described by the equilibrium particle and hole densities p{k) 
and p''(fc) of the charge degrees of freedom and the equilibrium string densities (t„(A) and cr^(A) of the spin degrees 
of freedom. Here n = 1,2,..., oo. The partition function Z = tr(e~^/"^) is defined by 

Z = W^(p,P^a„,a„^)e-^(''^'''^--")/^ (7) 



where the densities satisfy (U) with W{p, p''' , <Tn, cr'^) the number of states corresponding to the given densities. In- 
troducing the combinatorial entropy S = lnW{p, p'^ , cr^, a'^), the grand partition function is Z = e~'^^'^, where the 
Gibbs free energy G = E — pN — H{N^ — Ni)/2 — TS. Here p is the chemical potential and (N^) denotes the 
number of the particles with up (down) spin. Recall that = M. The energy per unit length is defined by 

/oo 
k^p{k)dk - m'^H. (8) 
-oo 

Here H is the external magnetic field and = {N-\ — Ni)/2 denotes the atomic magnetic momentum (where the 
Bohr magneton ps and the Lande factore are absorbed into the magnetic field H). The magnetization per unit length 
in the z-direction is thus given by 

-1 /"OO /"OO 

^'^n /^(^)* - II " / ^n(^)dA. (9) 

^ J — oo ^ J — oo 

Now the equilibrium states are determined by the minimization condition of the Gibbs free energy (30l . [s^ . i.e., 
the condition 5{E — pn — TS) — 0, which gives rise to the set of coupled nonlinear integral equations (the TBA 
equations) ^] 

eik) = - M - 2 g - ^ ^2 + (fc 1 fc,)2 HI + - r ^ a„(fc - A) . ln(l + r,-^ (A)) 



ln77„(A) - — + «n(A - k) * ln(l + e'^C^')/^) + ^ T„„(A - A') * ln(l + Vn'i^'))- (10) 

ri=l 

Here we have defined the dressed energy e{k) := T h\{p^ {k) / p{k)) with respect to the quasimomentum k. Similarly 
ry(A) := a'^{X)/a{X) with respect to the spin quasimomentum A. The dressed energy e(A;) plays the role of excitation 
energy measured from the Fermi energy level e(A;F) = 0, where fcp is the Fermi momentum. The pressure P{T, H) 
and free energy F(T, H) per unit length in the thermodynamic limit are given in terms of the dressed energy by 

rji pOO 

P{T,H) = — dk\n{l + e-<^^f^) (11) 

F{T,H) = pn / dfcln(l + e-'('=)/'^). (12) 

27r J^oo 
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At low temperatures and zero magnetic field, it is important to note that e ^^^^/'^ is negligibly small if the energy 
> 2/L(, where fi is the chemical potential. It follows that for strong coupling, i.e., for 7 ^ 1, the TBA equations 
([TP)) can be written as 

6(fc) = fc^ - /i - - ^'rf ;f ^ _ rf; /°° a4A) ln(l + ^'^(A)) 

1 f\\ niJ 47rP(T,g)a„(A) , nc.^ 

lnry„(A) = — + -^^ + 2^T„„(A-A)*ln(l + r;„ (A)). (13) 

?7 = 1 

In the above equations, we have only kept terms to order 1/c for the dressed enery e(fc). We have also made the variable 
changes A cA/2 and 77(A) ■q{c\/2). Thus the function a„(A) defined in ^ becomes a„(x) = ^ nz'^^i 7 which 
is Takahashi's notation for the ferromagnetic Heisenberg chain [sOl- It will be clearly seen from the TBA equations 
P3p that the low temperature behaviour of the strongly interacting two-component boson model is determined by 
the hard core Bose gas state and ferromagnetic spin wave fluctuations described by the ferromagnetic Heisenberg 
chain with an effective coupling strength J = 2P{T, H) jc > 0. We shall see through the TBA equations (fT3)) that 
the temperature dependent part of this effective coupling strength contributes to the free energy at 0{T^). At low 
temperatures J — > AE-p/S, where Ep is the Fermi energy. 

Following the analysis of the TBA equations for the ferromagnetic Heisenberg chain , the TBA equations ([13]) 
can be rewritten in terms of the function 771(A) in the form 

„ 2cP{T,H) 



where 



eik)^k'-^i + fxxx (T, H) (14) 

/oo 
dAs(A)ln(l+ 771(A)) (15) 
-00 

is the free energy of a ferromagnetic Heisenberg chain with coupling J w 2P{T, H) /c. Here l/s(A) = 4cosh(7rA/2). 
The function 771(A) is determined by the TBA equations for the ferromagnetic Heisenberg chain [30|, which are 

ln77i(A) = ^s(A)+s*ln(l + 772(A)) 

ln77„(A) = s*ln(l + 7;„_i(A))ln(l + 77„+i(A)). (16) 

For large n, the string solutions obey lim„^oo 7^~^ln77„(A) = H/T. 

So far we have separated the ferromagnetic spin state from the hard core gas phase via the TBA formalism. It is 
important to note that for the strong coupling regime the spin and charge degrees of freedom are coupled via the 
pressure. The spin velocity vanishes whereas the charge velocity tends to the Fermi velocity of noninteracting spinless 
fermions as 7 — > cxd. In this extreme case, the effective mass takes the maximum value m* — Nm, meaning that by 
moving one boson with down spin, one has to move all the particles with up spins [23.[25|. For finitely strong repulsion, 
a quadratic dispersion of spin wave excitations above the ferromagnetic ground state can be obtained from the Bethe 
ansatz solution ^ [24 . l25i] . In general spin charge separation is a typical phenomenon in interacting many-body 
systems [28j . 

Before moving on to discuss the solutions of the TBA we proceed to calculate the pressure (|lll) . Without loss 
of generality, we let e{k) = ^ - A{T, H), where A{T, H) ^ ^ + 2P{T, H)/c - fxxx{T, H). Hereafter we restore 
physical units in the thermodynamic properties. Integration by parts gives 



p(r, H) « r ^ttU, • (17) 

1 + e '^B-r 



2m 



The integral in ([TT]) can be calculated explicitly using Sommerfeld expansion [40] to give 



P(T,i/) « ^=^lA{T,H)i 

'jfTif. 3 
2 m 



KbT V 7tt^ f KbT 



A{T,H) J 640 \A{T,H) 



(18) 



Here the function A{T, H) contains P(T, H). From equation ([T5|) we can find a relation between the pressure P(T, H) 
and chemical potential /i by iteration. This provides a starting point to calculate the thermodynamics of the strongly 
interacting two-component Bose gas at low temperatures, including the zero temperature limit. 
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IV. THE GROUND STATE AND THE SPINLESS BOSE GAS 



Pure dynamical interaction drives the spinless Bose gas into distinct quantum phases of matter: from the quasi- 
Bose-Einstein condensate to the Tonks-Girardeau phase. As stated in the introduction, this elegantly simple Bethe 
ansatz solved ID quantum many-body system [9|, |10j is testable in experiments on trapped quantum gases of ultracold 
atoms [1, H B 0- -"-t is natural to expect that spin dynamics in the interacting two-component spinor Bose gases 
would lead to significantly different quantum effects than those of the spinless Bose gas. At T ^ 0, it is suitable to 
use the dressed energy formalism in the TBA equations (|16p . i.e., 

a (A) = 2ttJs{X) + s^^+{\) 

CnW = (19) 

where lim„^oo ^n/n = H. Here the dressed energy is defined as ^„(A) T\ni]n{X), with £,nW i^nW) denoting the 
dressed energy for ^„(A) > ('fn(A) < 0). The free energy is now 

fxxx{T,H) - -:,H + Y^ / a4A)C(A)dA 

= J In 2- / dXs{X)^tW- (20) 



For the ferromagnetic case J > 0, the solution of the TBA p9| is given by [s^l 

?n(A) = C(A) = 2^Ja„(A) + Hn (21) 

where n — l,2,...,oo. Here we see that for T — 0, fxxx = —11/2. Therefore it follows that the fully-polarized 
state forms a ferromagnetic ground state. In this way the TBA gives rise to a direct proof of the existence of the 
ferromagnetic ground state. 

In order to understand spin effects in the spinor Bose gas, we first discuss the low temperature behaviour of the 
spinless Bose gas. It was for this model that Yang and Yang introduced the TBA formalism, with result 

eik) . eOik) _ , - - -^^ ^.^(,-,O^ Ml + e-C^ )/-) (22) 



3 \ 7 / \ 37 

Po-^nV(l-^), Q^nn[l-'-]. (23) 



which is the special case of the TBA equation (|10|) for the spinor Bose gas. At zero temperature, one can obtain 
physical quantities, such as the ground state energy per unit length Eq, chemical potential ^, pressure Pq and the 
cut-off momentum Q. In the strong coupling limit the results for these quantities are 

1 3 2 / 4\ , 2 16 

Eq « —n TT 11 , /^o ~ '^"71' 1 — 



The macroscopic velocity is 

= W 9 _ , 

9?! m \ 7 

The low temperature thermodynamics can also be calculated directly from the pressure (|18p . For the spinless case, 
the function A{T, 0) « ^(1 -I- Ml^lMy Substituting A{T, 0) into equation (HH) and using the relation dP{T, 0)/9/i = n 
gives the chemical potential 



/ dPo n-KTi / 4\ 



Mo 



16\ f KbTV tt^ ( 32\ f KbTV 



37 j I MO ) ^ 36 57) ( MO j 



(25) 



Here /io ~ n^vr^ ^1 — coincides with the result given in equation In terms of the degenerate temperature 

T = KBT/Td the free energy per unit length follows from relation (fT2)l as 



F(r) « Eq 



1 -\ - —- (1 — 



(26) 
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The ground state energy Eq agrees with the result given in p3|) . The expression (|T8|) indeed provides a simple way to 
derive the thermodynamics. In addition the results (|^5|) and (pS)) obtained for the spinless Bose gas via the TBA are in 
good agreement with results derived from generalized exclusion statistics [4l|. The results (|25l) and (|26l) characterize 
the low temperature behaviour of the ID spinless Bose gas induced by the dynamical interaction in strong coupling 
regime. The specific heat Cy and the entropy S follow from the free energy ([26]) as 



C71 



TLd^F{T,0) _ NKbt 
TL dF{T, 0) 



NKbt^ 



6(1-^) 
NKbt 



dT 



6(1 



157r2(l- 
NKb 

12^ 

7 ' 



12) 

7 ' 



457r2(l 



(27) 
(28) 



which coincide with the results given in Refs. [4i|,|4j. Here the degenerate temperature r = KsT/Td, with Td 
The total energy per unit length follows from the relation E{T, 0) — F{T, 0) + ST, with result 



2m 



E{t) « Eo 



1 



T 

4^ 



1 + - 

7 



207r4 



1 + 1^ 

7 



The ground state energy Eq is as given in ([23|) . From ([26|) we see that 



F{T) = F(0) 



7rC(jrBr)' 



0(T2 



(29) 



(30) 



as expected from conformal field theory arguments for a critical system, i.e., for a system with massless excitations 
p3 | Here the central charge C = 1 and Vc is given by (I24p . Similarly the finite-size corrections [i^ are given by 



E{L,N)~Leoo = - 



hnCvc 
6L 



Oil/L^ 



(31) 



Here E{L, N) is the finite size ground state energy and Cqo is the energy per unit length in the thermodynamic limit. 

Furthermore, at low temperatures the strongly interacting spinless ID Bose gas can be viewed as a system of 
ideal particles obeying nonmutual generalized exclusion statistics (GES) with statistics parameter [ill, SE] a « 1 — 
2/7. These particles obey GES interpolating between bosons and fermions [H, [13, S^. For ID interacting many- 
body systems the pairwise dynamical interaction between identical particles is inextricably related to their statistical 
interaction through scattering. GES is thus the result of collective behaviour exhibited in ID quantum many body 
systems. From the GES approach, the free energy and the total energy per unit length are given by [4l| 



Fir) 
Fir) 



En 



Eo 



1 



T 

^2 



1 



7 

T' / 2 

47r2 ^ J 



27776 



1 



T 



607r4 
207r4 



1 



(32) 
(33) 



Here C(3) = E^^i V"'- 



V. LOW TEMPERATURE FERROMAGNETIC BEHAVIOUR 



The low temperature behaviour of the spinor Bose gas, triggered by the ferromagnetic spin-spin interaction, is 
intimately related to the thermodynamic behaviour of the ferromagnetic Heisenberg chain, which has been extensively 
studied via various methods, e.g., numerics [H, HOIj spin wave theory [3^ and the TBA approach with extrapolation 
[51I I . Although there has been a wide range of interest in the ferromagnetic Heisenberg chain, realization of ferromagnet 
chains are relatively rare [s^ . Most recent interest in the ferromagnetic Heisenberg chain has been from the perspective 
of string theory [53, [l^l • Obtaining exact analytical results for the thermodynamics of this model still provides a 
number of open challenges. For one, it is extremely hard to solve the infinitely many equations (|16p involved in the 
TBA. Nevertheless, Takahashi and his coworkers |36i, have given some results for the free energy and susceptibility 
which are generally accepted. Schlottmann [s^ has also predicted the leading order of the specific heat and zero-field 
susceptibility via analysis of the string solutions to the TBA equations p6|) . With the help of these known results for 
the ferromagnetic Heisenberg chain, we show here that, in the strong coupling regime, the ferromagnetic state induced 
by the internal spin-spin interaction significantly affects the low temperature behaviour of the two-component spinor 
Bose gas. 
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A. Paramagnet: 7 ^ oo 

We first consider tlie extreme case 7 00, or say 7 ^ 1/KbT. In tliis case tlie driving term in the TBA equations 
P6p vanishes as 7 00. Thus the string solutions are given by [soj 



sinh( 



{n+l)H 
2T 



sinh(^ 



(34) 



which are known as free spin solutions. In this case, the particles with down-spins are unable to exchange their 
positions with the particles with up-spins. The spins are thus frozen locally and the spin-spin exchange interaction 
vanishes. In this case, the statistical interaction is completely suppressed due to the strong repulsion. In this sense, 
for 7 — !■ 00 both the spinor Bose gas and the Fermi gas behave like a free spin paramagnet. In this case, the spins are 
very sensitive to external magnetic fields, with a small field able to polarize all atoms. From equation (|14p we obtain 



2 1,2 



e{k) 



2m 



i..Tln(2cosh^). 



(35) 



Under the condition 7 ^ 1/KbT and KbT < Td, the free energy gives 



F{T,h) 



3 2m' 



T 

4^ 



607r4 



3-r, , h 

In 2 cosh — 
2r 



(36) 



where we have set h — H/T^- 

It follows that in the strong coupling limit 7 — > cx) the magnetic properties of the two-component spinor Bose gas 
are those of su{2) free spins with a divergent susceptibility x ~ ak^t ~ ^SLnh^{H/{2KBT))) at low temperatures. 
The magnetization is = ^nta,nh{H/ {2KbT)). Figure [T] shows the free spin behaviour in the magnetization as a 
function of the magnetic field at different temperatures. 



The specific heat and entropy follow from equation (|36p as 



S 



\nKbt 
6 

\nKbt 
6 



NKbt^ NKbH 



2 r 



157r2 

NKbt^ 
457r2 



4t2 
NKb In 



1 - tanh^ ( — 
2r 



2 cosh — 
2t 



NKBh ^ ,fh 

tanh — 

2t \2t 



(37) 
(38) 



In the absence of external field {H = 0) the specific heat behaves like that of a free Fermi gas. However, the free spins 
make a contribution NKb In 2 to the free Fermi entropy. Figure [2] shows a plot of the specific heat in the presence of 
magnetic field. We see clearly here in the strong coupling limit that the specific heat is sensitive to the external field. 
We also note that our calculations differ from those in Refs [22|, [3^ where the spin degrees of freedom were ignored. 



B. Ferromagnetic: 1 7 < 1/KbT 

For finite temperatures the solutions of the TBA equations (fTH]) vary from the free spin solutions ([2T|) and ([34]) 
[30l . [3^ . The spin-spin exchange interactions are enhanced as the interaction strength 7 decreases from the strong 
coupling limit. In this section we explore the ferromagnetic behaviour of the spinor Bose gas with finitely strong 
interaction, i.e., in the regime 1 ^ 7 < 1/KbT, or more precisely speaking, 2p(t 0) ^ ^/KbT. In this regime, the 
known results for the free energy of the ferromagnetic Heisenberg chain at low temperatures obtained by Takahashi 
and colleagues are applicable to the ferromagnetic state associated with the free energy The result is (36l. Isij 



fxxxiT,0) « J 



-1.042 



KbT 
J 



KbT 
J 



0.9 



(39) 
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FIG. 1: (color online) Magnetization (normalized by the linear density n) vs magnetic field H (with Bohr magneton fiB ~ 1) 
for different temperatures (in units of Kb)- In the strong coupling limit 7 cx) two-component spinor bosons exhibit free spin 
behaviour at low temperatures. 



with effective coupling strength J w 2P{T,0)/c. At low temperature, the effective coupling approaches the constant 
value J ~ 4£'f/37 as a result of the temperature-dependent part in the coupling only making an 0((fcsT)^) contribu- 
tion to the free energy of the spinor Bose gas. Schlottmann [35| calculated the leading terms of the specific heat and 
free energy for the ferromagnetic Heisenberg chain by approximating the infinite set of TBA equations Using 
Schlottmann's method we find from the TBA equations p6p that the leading order in the free energy is proportional 
to (KbT)^ . With the help of Takahashi's result for the ferromagnetic Heisenberg chain (|39p . we may calculate the 
thermodynamics of the spinor Bose gas within finitely strong interaction, which we now do. 
Substituting A{T, 0) = + 2P(T, 0)/c - fxxx{T, 0) into equation dH]) we have 
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where jl = ^ — ,fxxx{T,0). For the regime 2p{t 0) ^ ^/^bT calculation of the chemical potential via the relation 
aP(r,0)/9/i = n gives 
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where 
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Some algebra gives the free energy F(t, 0) and total energy E{t, 0) per unit length 
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FIG. 2: (color online) Specific heat normalized by NKb vs degenerate temperature t = KbT /Td- The straight dashed line is 
the specific heat (|27[) for the spinless Bose gas with 7 = 20. The solid line is the specific heat for the free fermions. The curves 
show the specific heat (I37[) of the spinor Bose gas with 7 cx) for finite field term h. 
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We see explicitly that the ferromagnetic behaviour of the spin exchange interaction dominates the thermodynamics 
of the strongly interacting two-component Bose gas. 

It is of interest to note that the leading temperature-dependent term in the energy expressions is 0{T^^'^) for the 
spinor Bose gas compared to 0{T'^) for the spinless Bose gas (|26l) . In the strong interaction regime 1 <C 7 < ^iq/KbT 
this leads to quantitatively different low temperature behaviour compared to the spinless Bose gas. Figure [3] shows 
the total energy and the free energy in the strong interaction regime as a function of temperature for the spinor Bose 
gas, the spinless Bose gas and an ideal gas obeying GES. Significantly different characteristics of low temperature 
behaviour for the spinor Bose gas and the spinless Bose gas are depicted. For the spinless Bose gas with strong 
coupling, the low temperature thermodynamics is known to coincide with that of ideal particles obeying nonmutual 
GES with statistics parameter a « 1 — 2/7. Figure [3] hints that the strongly interacting spinor Bose gas might also 
be equivalent to a gas of ideal particles obeying nonmutual GES. It is an open question as to what the statistics 
parameter a for the spinor Bose gas might be. 



To conclude this section the specific heat and entropy are given by 
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FIG. 3: (color online) Free energy F{T,0) and total energy E{T,0) per unit length (normalized by nE-p = ■^T^^n'^) as a 
function of the degenerate temperature r — KbT /Td for different coupling strengths for both the two-component spinor and 
spinless ID Bose gases. The thick solid and dashed lines are the energies (|42p and (|43|) of the spinor Bose gas. The flatter thin 
solid and dashed lines are the results (|26|l and (|29p for the spinless Bose gas. For further comparison the symbols show the 
energies (|32p and (|33p derived for the spinless Bose gas from an ideal gas of particles obeying nonmutual GES In general 
the curves reveal the universal characteristics of the energies at low temperatures. 



S 1.042 X 3\/3(7t)5 3(7r) 0.9 x 15\/3(7r)i t 

^ 4(1 -fK 2(l-f)7r2+ 16(1- +6(l-f) ^ ' 

which differ significantly from the corresponding spinless Bose gas results (|27|) and ([28]) and also from free spin case 
(157)) and ([55]) . Figure m shows the specific heat as a function of the temperature. The specific heat exponent c« ~ T~° 
indicates that a = —0.5 for spinor Bose gas for the regime 1 ^ 7 < iiq/KbT whereas a = — 1 for the Lieb-Lininger 
gas. For 7 ^ 1/ KbT, the spinor Bose gas and the Lieb-Lininger gas both have a = —1 for the absence of the external 
field. 



C. Susceptibility: K 7 < 1/KbT 

Now we consider the effect of a small external field {H ^ T) within the regime 1 <C 7 < 1/ KbT. Here we adapt 
the known free energy result for the ferromagnetic chain in the presence of an external field [36l . [5l| , namely 
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with fxxx{T, 0) as given in ([39]) . Repeating the procedure of the previous section with the free energy ((46|l gives the 
result 
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with F(r,0) as given in ((i^ . 
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FIG. 4: (color online) Comparison between the specific heat (normalized by NKb) of the two-component spinor Bose gas and 
the spinless Bose gas as a function of the temperature parameter 7r at difTerent coupling strengths. The upper set of curves 
are obtained from the spinor Bose gas specific heat (|44|) . The lower set of curves follow from the spinless Bose gas result p7|) . 
In general the curves highlight the distinction between the temperature-dependent behaviour of the spinor and spinless Bose 
gases. 



For small magnetic field (h < r) the susceptibility per unit length 
d^F{T,H 
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(48) 



is indeed greater than that of free spins, for which x ~ n/{4:KBT). This result is also consistent with Eisenberg and 
Lieb's general argument for polarized spinor bosons psj . Figure [HI shows the zero-field susceptibility as a function of 
temperature for different values of the interaction strength. The susceptibility diverges as r — > 0. The susceptibility 
decreases with increasing interaction towards the free spin paramagnetic susceptibility is the lowest curve (solid line) . 
The susceptibility exponent defined by x T^'' is b = 2 for the regime 7 < ^jlo/KbT with 6 = 1 for the paramagnet. 



VI. LOCAL PAIR CORRELATION 



The local pair correlation function for the ID Bose gas has been determined experimentally in a gas * Rb atoms 
as a function of the interaction strength by measuring photoassociation rates [55|. In general local two-particle 
correlations can be used to study phase coherence behaviour and classify various finite temperature regimes in ID 
interacting quantum gases 56, 57| . In the grand canonical description, two-particle pair correlation are given in terms 
of the field operator ^ and the free energy f{j,T) by [2^, [sl] 

gi^H0):^i^^^^^^,) = ^(^Ih2ll) . (49) 

At zero temperature the local pair correlation is g^"^^ « 1 for the weakly interacting Bose gas. On the other hand, 
g^"^") — > as 7 increases into the Tonks-Girardeau regime, indicative of free fermionic behaviour. In this regime, the long 




temperature x 



FIG. 5: (color online) Susceptibility (normalized by n/Td) as a function of the degenerate temperature r for different values of 
interaction strength 7. In the strong coupling limit 7 — > 00 the two-component spinor bosons behave like a paramagnet with 
susceptibility x ^ T^^ . For finitely strong interaction, the ferromagnetic susceptibility behaves as x ~ . 



range behaviour is characterized by the one-body correlation function g^^\x) = (4'^(a;)^'(0)) oc corresponding 
to the momentum distribution n{p) cx l/^/p. Here p is the momentum. More generally, in terms of the Luttinger 
parameter K 58], g'^^\x) oc l/x^^^^ and n{p) cx In the weak coupling limit K w 7r/^/7, which leads to 

a power-law decay in the one point correlation. 

The pair correlation function for the spinless Bose gas 

9^'\0)^'^(l + ^,+^\ (50) 
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follows from the TBA result (|26p. which coincides with the result given in Refs. [5^,[53|. It is evident that the dynamical 
interaction dramatically reduces pair correlation due to decoherence between individual wave functions of colliding 
particles. On the other hand, increasing temperature slowly enhances local pair correlation. At temperatures r <C 1 
(T <g; Td), the local pair correlation approaches free Fermi behaviour as 7 — > 00, as was quantitatively demonstrated 
in the experimental observations of the pair correlation function for a gas of interacting *^Rb atoms confined to ID 

For the two-component spinor Bose gas considered here, the ferromagnetic spin-spin exchange interaction results 
in a different temperature-dependent pair correlation function. For this model, the local pair correlation function 



47r2 
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(51) 



follows from equation (|42| in the regime ^q/ KbT > 7 3> 1. We see for the spinor Bose gas the local pair correlation 
again quickly decays with respect to the dynamical interaction 7. Figure ^ shows the local pair correlation for the 
spinless and the spinor Bose gases as a function of the interaction strength at different temperatures. In contrast 
to the spinless Bose gas, where the temperature enhances local pair correlation due to thermal fluctuations, for the 




FIG. 6: (color online) Local pair correlation g^^^ as a function of interaction strength 7 at different temperatures. For the 
spinless Bose gas, the pair correlation function increases with increasing temperature. In contrast, for the two-component spinor 
Bose gas the pair correlation function decreases with increasing temperature as a result of the temperature induced spin-spin 
ferromagnetic exchange interaction. 



spinor Bose gas the local pair correlation decreases with increasing temperature due to the ferromagnetic spin-spin 
exchange interaction. 



VII. CONCLUSION 



In this paper we have studied the thermodynamics of the integrable ID two-component Bose gas via the thermo- 
dynamic Bethe ansatz. Analytic low temperature results were obtained for the free energy, total energy, specific heat, 
entropy, pressure, susceptibility and pair correlation function in the strongly interacting regimes 7 3> 1/KbT and 
1 «C 7 < ^Iq/KbT. Where appropriate, comparison was made with corresponding thermodynamic properties of the 
integrable ID spinless Bose gas. Our key finding is that the temperature induced ferromagnetic spin-spin exchange 
interaction triggers a number of novel quantum effects in the thermodynamic properties of the spinor Bose gas at 
low temperatures. In the regime 1 <C 7 < ^lq/KbT the specific heat exponent for the spinor Bose gas following from 
Cv ~ T^/^ is different to that of the spinless Bose gas for which Cv ^ T. In this regime, the susceptibility exponent 
is given by x '--^ -^i^^q]^ exceeds that of free spin paramagnet for which x ^ T~^. In contrast to the spinless 
Bose gas, where the pair correlation function increases with increasing temperature, the two-component spinor Bose 
gas pair correlation function decreases with increasing temperature as a result of the temperature induced spin-spin 
ferromagnetic exchange interaction. 

In general these exact results should be relevant to understanding ferromagnetic behaviour and spin effects in two- 
component spinor Bose gases of cold atoms, for which the interaction strength can in principle be tuned. However, the 
introduction of precise thermometry into these systems to measure universal temperature dependent effects provides 
a number of challenges. However, this is a worthwhile goal. As we have seen, for strong coupling the characteristics of 
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the thermodynamics of the spinor Bose gas at low temperatures are described by an effective ferromagnetic Heisenberg 
spin chain. There is a remarkable three-way correspondence [ssl [s^ between the ferromagnetic Heisenberg chain, 
the limit of weakly coupled planar A/" = 4 supersymmetric Yang-Mills theory and the limit of free strings on AdSs x 
S^. This triality between guage theory, string theory and the thermodynamics of the ferromagnetic Heisenberg 
chain has recently been used to calculate the Hagcdorn temperature of the string theory in agreement with the 
Hagedorn/deconfinement temperature calculated on the guage theory side [54|- It now appears that we can add a 
further connection with the thermodynamics of the strongly interacting two component spinor Bose gas. 
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